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Kink Structure in the Quasiparticle Band of Doped Hubbard Systems 
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By making use of the self-consistent projection operator method with high-momentum and 
high-energy resolutions, we find a kink structure in the quasiparticle excitation spectrum of the 
two-dimensional Hubbard model in the underdoped regime. The kink is caused by a mixing 
between the quasiparticle state and excitations with short-range antiferromagnetic order. We 
suggest that this might be the origin of the strong concentration dependence of the 'kink' found 
in Laa-^Sr^CuCU (a: = 0.03 - 0.07). 

KEYWORDS: kink, quasiparticle excitations, Hubbard model, ARPES, cuprate, LSCO 



Recent high-resolution photoemission experiments 
show that there is a well-defined 'fcmfc' in the quasipar- 
ticle band dispersion of high-T c cuprates, whose energy 
scale is Wkink = 60 — 70 meV in both the normal and 
superconducting states. 1-3 The kink along the nodal di- 
rection was found to have a universal feature, 2 i.e., the 
Fermi velocity vf in the low-energy region (\u>\ < Wkink) 
is not sensitive to the type of cuprates, doping concen- 
tration, and isotope substitution, 4 although vf in the 
high-energy regime (|o;| > i^kink) strongly depends on 
the latter. Various theoretical explanations for the kink 
have been attempted from two different points of view. 
One type of theory relies on the coupling of an electronic 
quasiparticle to the spin fluctuation resonance mode ob- 
served in inelastic neutron scattering experiments. 5, 6 An- 
other relies on a coupling to phonons, particularly to the 
longitudinal optical phonon mode found in neutron ex- 
periments. 2 ' 3 The latter approach raised again the fun- 
damental question on the mechanism underlying high-T c 
superconductivity in cuprates, i.e., electron- or phonon- 
mediated. Although the reliability of these theories is un- 
der debate and the improvements of these theories are in 
progress, 7,8 it has not yet been seriously studied whether 
or not the kink in the quasiparticle state can be solely 
of electronic origin. We deal with this problem in the 
present letter and report the appearance of a kink due to 
long-range electron correlations in the two-dimensional 
(2D) Hubbard model for small doping concentrations. 

The difficulty in the present problem is that the per- 
turbation approach is not applicable to cuprates be- 
cause of strong electron correlations. Therefore, ad- 
vanced theories such as the Lanczos method, 9 the quan- 
tum Monte-Carlo (QMC) method, 10-12 and dynamical 
cluster approximation (DCA), 13,14 have been applied to 
the cuprate system. They clarified the global structure 
of the single-particle excitation spectrum in the 2D Hub- 
bard model. A detailed structure of the low-energy ex- 
citations at low temperatures, however, has not been 
derived because of the limited resolutions in both mo- 
mentum and energy and the limited range of intersite 
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electron correlations inherent in the cluster approaches. 
In particular, a high resolution is indispensable for theo- 
retical investigations of the kink. In order to overcome 
these difficulties, we have recently developed the self- 
consistent projection operator method (SCPM) on the 
basis of the retarded Green function. 15, 16 The method 
is suitable for the present purpose because it allows for 
calculating self-energy with high-momentum and high- 
energy resolutions, and self-consistently taking into ac- 
count the long-range intersite correlations by using an 
incremental cluster expansion and an effective medium. 
The results presented here were obtained by this method. 

We apply the 2D Hubbard model on the square lat- 
tice using the nearest-neighbor electron hopping param- 
eter t and intra-atomic Coulomb interaction parameter 
U. The single-particle excitation energy spectrum is ob- 
tained from the retarded Green function 

Here, z — to + i5, where 5 is a positive infinitesimal 
number, and e k is the Hartree-Fock one-electron energy 
dispersion measured from Fermi energy. In the SCPM, 
the momentum-dependent self-energy Ai.(z) is calcu- 
lated from the Fourier transform of nonlocal memory 
functions My as 



A k (z) = U 2 ^2 M M Z ) exp(ifc • Rj) 



(2) 



Note that Rj is the position vector of site j. High- 
momentum and high-energy resolutions are achieved by 
taking into account the off-diagonal terms Mij(z) up to 
infinity. 

We calculate Mij(z) by means of an incremental 
cluster expansion in an effective medium with a co- 
herent potential Within the two-site approxima- 
tion, the Mij(z) are given by Mu(z) = m\1\z) + 
E^i(M^(z)-M^(z)) and M^{z) = M^{z). 
M^f (z) and Mj^j (z) are the matrix elements of the clus- 
ter memory matrices defined by M$(z) 
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Fig. 1. Excitation spectra along high symmetry line calculated at 
hole concentration <5 h = 0.14, U = 8, and T = 0. T = (0, 0), X = 
(0, it), and M = (it, it) in units of the lattice constant. The energy 
unit is chosen so that the nearest-neighbor transfer integral is 
one. Open circles with error bars are the results obtained by 
the QMC method 12 at T = 0.33. The dashed curve shows the 
Hartrce-Fock contribution e^. 



(for c = i) or (2x2) (c = (ij)) cluster Liouvillcan, 
whose diagonal matrix elements are given by L^(z) = 
J7(l-2(n i _ CT ))/[(n i _ CT )(l-(nj_ -))] using an average elec- 
tron number {n^) with spin a on site i. The screened 
memory matrix element M-^(z) is obtained using renor- 
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Fig. 2. Contour map of excitation spectrum along nodal direction 
at <5 h = 0.01 



malized perturbation theory 15 as 



z - e - e' + e" 



MV{z) = A, 

! . — I — ! -t- : 

(3) 

where A u = [(n i _ (r )(l - (n i _ (r })]/[(n i _ (7 } c (l - (n;_ CT ) c )] 
and A^j = 1. Furthermore, (rii a ) c — J depf^ (e)/(e) 
with the Fermi distribution function /(e). Moreover, 
X(e,e',e") c = f(-e)f(-e')f(e") + f( e )f(e')f(-e")- The 
matrix pf^ (e) is the density of states of a system with 
an empty site i (or sites i and j) embedded in a medium 
E(z). The latter is determined self-consistently by means 
of the coherent potential approximation (CPA) equa- 
tion. 18 We note that the self-energy obtained using eq. 
(3) yields the second-order perturbation theory in the 
limit of a small U and reduces to the exact result in the 
limit of a large U. In the intermediate coupling regime, 
we find a quantitative agreement of quasi-particle bands 
between the SCPM and QMC method as will be shown 
below. 

In the numerical calculations, we adopted the two- 
site approximation and took into account intersite cor- 
relations up to the 50th nearest neighbor site. More- 
over, we assumed excitations from the paramagnetic 
ground state (i.e., the normal state). Figure 1 shows the 
momentum-dependent excitation spectrum of a doped 
Hubbard system for an intermediate Coulomb interac- 
tion strength J7 = 8|t|atT = and a hole concen- 
tration <5h = 0.14. The excitations consist of an upper 
Hubbard band around the M (7r,7r) point, incoherent 
excitations at the T point with an energy u> ~ A\t\ as 
traces of the lower Hubbard band, and the quasi-particle 
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Fig. 3. Contour map of excitation spectrum along nodal direction 
at <5 h = 0.05 



band near the Fermi level. We emphasize that the quasi- 
particle band agrees well with that obtained by the QMC 
method 12 for the same value of U. We do not find, how- 
ever, any kink in the quasi-particle dispersion in the opti- 
mal doping regime, although the resolution of our results 
is much improved as compared with that obtained by the 
QMC method. 

We have examined the details of momentum depen- 
dent excitations along the nodal direction ((0, 0) — (ir, it)) 
from the underdoped to the overdoped regime. Figure 2 
shows the results for a doping concentration of 5^ = 0.01. 
We find here a kink at |fe| = 0.57T and Wkink = — 0.8|t|. 
The kink becomes weaker with increasing doping concen- 
tration and the velocity ratio of the high-energy part to 
the low-energy one (v' f /vf) at the kink position becomes 
smaller (e.g., v' F /vp — 1.8 for Sh — 0.01, and v' f /vf = 1.5 
for Sh — 0.02), while the position hardly changes with 5^. 
As shown in Fig. 3, the kink disappears for 5^ = 0.05. 
This occurs together with a collapse of the lower Hub- 
bard band. 

It should be noted that in Fig. 2 the flat band exists 
at w = ±2.0|t| because of the excitations due to short- 



J. Phys. Soc. Jpn. 



Letter 



Author Name 3 



2.5 



1.5( 



0.5 



oo 



OoOOOOOqOO 



oo 



oo 



(0,0)-(7t,7t) U=8 



0.05 0.1 

5h 



0.15 



0.2 



Fig. 4. Fermi velocity along nodal direction as function of hole 
concentration. Closed circles denote the velocity in the marginal 
Fermi liquid state, while open circles indicate the velocity in the 
normal Fermi liquid state. The dashed line is the result for a 
nonintcracting system. 



range magnetic order (i.e., a precursor of the gap for- 
mation due to antiferromagnetic correlations). 11 ' 13 The 
kink at Wkink = — 0.8|t| is caused by a mixing between 
the quasiparticle state and magnetic excitations. Re- 
cently, we have reported that a marginal Fermi liquid 17 
(MFL)-like behavior is found in the undcrdoped region 
[Sh ^5 0.03) away from half-filling because of a pinning 
of Fermi energy to the van Hove anomaly due to a trans- 
fer of spectral weight from the lower Hubbard band to 
the upper one. 18 Antiferromagnetic correlations should 
be enhanced in this region because of nesting. The kink 
behavior appears in this region 5h < 0.03. On the other 
hand, the lower Hubbard band collapses when 5^ ^ 0.03. 
The MFL-like state with antiferromagnetic correlations 
changes to a normal Fermi liquid state and the kink dis- 
appears, as noted in Fig. 3. The calculated Fermi velocity 
along the T-M line is presented in Fig. 4. The velocity 
shows a weak concentration dependence, and the renor- 
malization factor of the Fermi velocity is about 1.8 from 
the underdoped region to the overdoped one. 

We have also examined in detail excitations along T- 
X-M. Because of the van Hove singularity, the quasipar- 
ticle band around the X point is quite flat. We do not 
find a kink behavior in this flat band region. The mix- 
ing between the quasiparticle band and the band of the 
magnetic excitations takes place away from the linear 
dispersion regime near the Fermi level, i.e., (\k\,uj) = 
(0.57T, —1.2) on the T-X line. Because both bands are 
flat, it is not clear whether this region contains a kink. 

The present model is too simple for attempting a 
comparison with the experimental data. Nevertheless, 
it is plausible that the kink is enhanced with decreas- 
ing doping concentration because of the development 
of short-range antiferromagnetic order. We speculate 
that the strong concentration dependence of the kink 
in L^-xSr^CuCU (LSCO) 3 ' 19 may be caused by the 



present mechanism. In fact, we obtain for the charac- 
teristic kink energy Wkink = 70 mcV when we choose 
the transfer integral so that the calculated Fermi veloc- 
ity along the nodal direction agrees with the observed 
one (1.8 eV-A). The value for Wkink agrees well with the 
experimental one (60-70 meV). 19 

In summary, we have investigated the quasiparticle 
band in the doped 2D Hubbard model on the basis 
of the SCPM with high-momentum and high-energy 
resolutions. We find a kink along the nodal direction 
(0,0) — (tt, 7r) in the range of doping concentrations 
< Sh ^ 0.03 where the MFL behavior persists. It is 
caused by a mixing between the quasiparticle excitations 
and magnetic excitations with short-range antiferromag- 
netic order. The kink decays rapidly with the decrease 
in antiferromagnetic correlations. We speculate that the 
kink in the undcrdoped regime of LSCO may be due to 
the present mechanism. 
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